Spin Susceptibility of Noncentrosymmetric Heavy-fermion Superconductor CeIrSi3 
under Pressure: 29 Si-Knight Shift Study on Single Crystal 
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We report 29 Si-NMR study on a single crystal of the heavy-fermion superconductor CeIrSi3 with- 
out an inversion symmetry along the c-axis. The 29 Si-Knight shift measurements under pressure 
have revealed that the spin susceptibility for the afe-plane decreases slightly below T c , whereas along 
the c-axis it does not change at all. The result can be accounted for by the spin susceptibility in 
the superconducting state being dominated by the strong antisymmetric (Rashba-type) spin-orbit 
interaction that originates from the absence of an inversion center along the c-axis and it being much 
larger than superconducting condensation energy. This is the first observation which exhibits an 
anisotropy of the spin susceptibility below T c in the noncentrosymmetric superconductor dominated 
by strong Rashba-type spin-orbit interaction. 
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Recent discoveries of superconductors that lack inver- 
sion symmetry have aroused great j^tCTejst J)ecause tliey 
uite unique prop ert iestjj, 

II 

the relationship between spatial symmetry and the spin 
state of Cooper-pairs is expected to be broken in those 
superconductors, which would make the parity of super- 
conducting state mixed between even and odd parities 
due to the antisymmetric spin-orbit coupling (ASOC), 
thereby leading to a two-component order parameter 
composed of spin-singlet and spin-triplet Cooper pair- 
ing states [H, H, 01, 0, H, IE 0, [§]• The heavy fermion 
(HF) superconductivity (SC) without inversion symme- 
try was first discovered in CePtsSi [9(. Subsequently, 
pressure (P)-induced SC was also observed in Ce-based 
HF antiferromagnets CeRhSi3 0, CeIrSi 3 jll], and 
CeCoGe3 jl2j , in which the inversion symmetry along 
the c-axis is broken. It is noteworthy that the SC tran- 
sition temperature T c = 1.6 K for CeIrSi3 under P= 2.6 
GPa is relatively high among the Ce-based HF supercon- 
ductors, which may be attributed to the strong coupling 
effect |13[ derived from the presence of strong antifer- 
romagnetic (AFM) spin fluctuations [3]. The most re- 
markable feature in the SC state of CeIrSi3 and CeRhSi3 
is that the upper critical field (i? C 2) exceeds 30 T along 
the c-axis, which is much higher than a Pauli limiting 
field III, 1(| 12 1- This kind of behavior has never been 
observed in other spin-singlet HF superconductors, such 
as CeCoIns, Celrlns, CeCu2Si2, and so on. Actually, 
the ASOC of CeIrSi 3 has been estimated to be 1000 K 
from the Fermi surface splitting in LaIrSi3 [15fl , which 
is much larger than a superconducting condensation en- 
ergy scaled to fee^c ~ 1.6 K. In this context, a novel SC 
state in noncentrosymmetric crystal structure would be 
expected for CcIrSi3 and CeRhSi3. In general, the mea- 
surement of Knight shift by means of Nuclear Magnetic 
Resonance (NMR) can provide a clue to determine the 



spin state of the Cooper pair [21 1. 

In this Letter, we report on the 29 Si-Knight shift mea- 
surement in the SC state of CeIrSi3 under P=2.8 GPa. 
Its spin component slightly decreases below T c for the 
a6-plane, whereas it does not change along the c-axis. 
This result reveals that the spin susceptibility Xs(T) for 
noncentrosymmetric superconductors without inversion 
symmetry is dominated by the anisotropy of the van- 
Vleck-like spin susceptibility originating from the tran- 
sition between the Fermi surfaces split by the ASOC, 
which makes it difficult to identify the admixture of the 
Cooper-pairing state between even and odd parity. 

A single crystal of CeirSi3 enriched by 29 Si isotope was 
grown using a tetra-arc furnace, as described in the pre- 
vious literatures 
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Il5l |. A small piece of the single 
crystal, size of which is 1.3xl.3x5mm, was cut out from 
the ingot of a large single crystal by using the x-ray Raue 
method. Hydrostatic pressure was applied by utilizing a 
NiCrAl-BeCu piston-cylinder cell filled with an Si-based 
organic liquid as the pressure-transmitting medium. To 
calibrate the pressure at low temperatures, the shift in 
T c of Sn metal under P was monitored by its resistiv- 
ity. 29 Si-NMR measurements have been made at a mag- 
netic field H =1.326 T. In the experiment of SC state 
under pressure, the pressure cell mounted on the mixing 
chamber was rotated within the plane including [110] and 
[001] directions, and fixed to the [110] or [001] directions 
by measuring the angle dependence of 29 Si-NMR spec- 
tra. The single crystal was adjusted with an angle of less 
than ±4° between the direction of H and [110] or [001] 
directions, as described later. 

FigureQJb) shows 29 Si-NMR spectra of the single crys- 
tal under P = 2.8 GPa at T = 1.1 K and H =1.326 T 
parallel to [110] and [001] directions. The lincwidth of 
each spectrum is as narrow as ~5 kHz, which is much 
narrower than that for the polycrystal shown in Fig.QJa), 
ensuring the high quality of the single crystal. There are 
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FIG. 1: (Color online) 29 Si-NMR spectra at H = 1.326 T and T = 
1.1 K for (a) the polycrystal and (b) the single crystal of CeIrSi3 
under P = 2.8 GPa. The NMR spectra at two inequivalent Si sites 
(denoted as Si(l) and Si(2)) are well separated in the single crystal. 

two crystallographically inequivalent Si sites denoted as 
Si(l) and Si(2), as illustrated in the inset of Fig.QJb). I n 
H || [110], the NMR spectra arising from Si(l) and Si(2) 
are well separated, as seen in the upper part of Fig.QJb). 
Since the number of Si(l) is twice as much as that of 
Si(2), the spectrum with twin peaks at lower frequencies 
is assigned to Si(l) and hence the spectrum with a sin- 
gle peak to Si(2). The twin peaks at Si(l) arise because 
H is not exactly set along the [110] direction, but with 
approximately ±4° from it. Note that the hyperfine field 
at Si(l), A$° ] (l) and hence Knight shift, K(l)^ is 
highly anisotropic within the a6-plane, originating from 
the pseudo-dipole field from 3p orbitals. By contrast, 
the single peak at Si(2) reveals that an anisotropy in the 
Knight shift X(2)t 110 l within the a6-plane is much smaller 
than that in ^(l)! 110 !. When H is applied to the [001] 
direction, the spectra at Si(l) and Si(2) overlap because 
the Knight shift (id 001 !) at both sites is almost the same. 
In general, Knight shift is given by, 

K(T)=K s (T) + K OTh + K dia 

where K S (T) and K m \^ are the spin and orbital part of 
Knight shift, and iidia is the SC diamagnetic shift. Note 
that K S (T) = j4hfXs(r) depends on T, where Am and 
Xs(T) are hyperfine coupling constant and spin suscep- 
tibility. K OI b is related to the usual van Vleck suscepti- 
bilities and hence does not depend on T. Figure [2Ja) 
shows the temperature (T) dependencies of K(l)^ 110 \ 
if(2)[ n °] and in the normal state up to 170 K 

under P = 2.8 GPa. The K(2)^ increases upon cool- 
ing, which is scaled to the bulk susceptibility x{T) at 
ambient pressure (P = 0) [15]. In contrast, i^(l)[ 110 l 
and if! 001 ] decrease upon cooling because Ams are neg- 
ative. Note that if(l)I U0 ], which was estimated from an 



FIG. 2: (Color online) (a) T dependence of the 29 Si-Knight shift 
in the normal state under P = 2.8 GPa. (b) Plots of if(l)[ 110 l and 
j^[00l] versus iC(2)f 110 l as the implicit parameter of T. Slopes in 
this plot give rise to the ratios of the hyperfine coupling constants, 
^li 101 WML" ' ( 2 ) - -0-45 and (X)/A^ (2) « -0.16. 



average value of shifts at the twin peaks in the spec- 
tra of Si(l) in Fig. [ljb), is in accord with the shift 
at the field along the [110] direction. Because of the 
lack of x(T) data under P, we could not make a plot 
of K versus x to deduce Am and K^. However, as 
shown in Fig. [H^b), respective plots of X(l)[ 110 ] and 
Jf[° 01 ] versus K{2)^ lltfi as the implicit parameter of T 
give rise to the ratios of A hf (1) [110] Mm (2) [no1 « -0.45 
and A M ° 1] /Am(2) [110] » -0.16, from their slopes of ten- 
tative linear lines that were indicated in the figure. 

Next, we deal with the Knight shift measurements in 
the SC state. Figures [3]Ja-d) indicate the T dependencies 
of 29 Si-NMR spectra and the Knight shifts below 1.5 
K under P = 2.8 GPa. The T C {H) = 1.2 K in the 
present experiment at P —2.8 GPa was determined 
from the result that the inhomogeneity of the field in 
the SC mixed state increases the NMR spectral width 
below 1.2 K as presented in the inset of Fig. [3jd). In 
H || [110], J FC(2)[ 11 °] decreases below T c , pointing to 
the decrease of Xs(T). Owing to the opposite sign 
of Am, i4T(l)[ 110 l slightly increases. Using the rela- 
tions of X(2)[ 110 l = ^ hf (2)[ 110 lxs(T)+^ orb (2)+X dia 

and K(l)^ = AM(l) [110] Xs(T)+K mh (l)+K d - m 
where -Kdia is identical at Si(l) and Si(2); and 
A h f(l)[ 11 °]/^hf(2) [1101 w -0.45, K dia was evaluated 
to be ~0.006% {H dla ~0.78 Oe in the field scale). Since 
an estimation of if rb at P = 2.8 GPa is difficult due to 
that a K-x plot is not available, ^ or b(2) = at P = 2.8 
GPa is assumed as well as at P = [22|]. A broken 
curve in Fig. [3](c) is a calculated curve for a case where 
Xs(T) would be reduced below T C (H) if a SC had a 
line- node gap with 2A/ksT c ~ 6 and the residual DOS 
of 37% at the Fermi level(.E'p), which was experimentally 
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FIG. 3: (Color online) T dependence of the 29 Si-NMR, spectra at 
H = 1.326 T parallel to (a) [110] and (b) [001] directions below 
1.3 K. The T dependencies of JC(l)[ 110 i and A'(2)[ 110 ' along (c) 
[110] and id 001 ! along (d) [001] directions. The inset of (d) is 
the T dependence of full-width at half-maximum (FWHM) of the 
spectrum for H || [001], which enables us to determine T C (H) = 1.2 
K. Broken curves are calculations for a case of Xs(T) — > at low T 
limit, assuming the residual DOS of 37% at the Ep deduced from 
the 1/Ti measurement on the polycrystal [13 . 



deduced from the 1/Ti measurement on the polycrys- 
tal |14| |. Likewise, by assuming the same model and 
by using the ratios of A M (l) [110] Mhf (2)[ 110 1 « -0.45 
and A^^'/Ahf (2)[ 110 1 w —0.16 obtained experimentally 
the increasing events below T c would be expected in 
-ft's(l)' 110 ' and if] 001 ', as shown by the broken curves 
in the figure. In H || [001], if! 001 ! should be increased 
if Xs(T) [001] were reduced below T C (H), since A [ °° 1] 
is negative, as indicated in Fig. [3^d). However, as 
decreasing T below T C (JT) = 1.2 K, slightly 
decreases by 0.006%, which is comparable to along 
the [110] direction. If we assume that i^dia for H || [001] 
is nearly the same as for H || [110], it might be expected 
that if] 001 ' does not change below T c . 

Figure H shows the plots of K S (T)/K S (T C ) versus 
T/T c for H || [110] and [001]. It is remarkable 
that (K S (T) I K S (T C ))[ 110 ] slightly decreases, whereas 
(K s (T)/K s (T c ))l° 0X i does not decrease at all. Note that 
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FIG. 4: (Color online) The T dependences of the plots of 
K s (T)/K a (T c ) versus T/T c for H || [110] and [001] in the SC state. 
A broken curve is a calculation for a case of Xs{T) — > at low T 
limit, assuming the residual DOS of 37% at the E F [3(See text). 



when usual superconductors with an inversion center are 
in a spin-singlet regime, Xs(T) decreases to zero gener- 
ally regardless of crystal directions, and when those are 
in a spin-triplet regime, Xs(T) for H _L d — vector stays 
constant, but that for H\\d — vector decreases to zero. 
It is known that the K S (T) in some spin-singlet SC re- 
mains finite at low T limit in association with an impu- 
rity and/or defects existing inevitably [3]. As a matter 
of fact, the present result cannot be reproduced by the 
simulation for a case of Xs(T) — > at low T limit, even 
though taking into account the residual DOS of 37% at 
the Ep due to the impurity effect [14|. as shown by a bro- 
ken curve in Fig. [H In this context, it is anticipated that 
the present results for CcIrSi3 with no inversion symme- 
try differ from these behaviours in the centrosymmetric 
superconductors. 

In inversion-symmetry-broken systems, it has been 
theoretically shown that Xs (T) consists of the Pauli term 



Pauli/ 



X'""*'(T) and the van-Vleck-like term x VV (T) [7]; The 
former originates from the Pauli paramagnetism due to 
the spin distribution on the intrabands. On the other 
hand, the latter originates from the transition between 
the Fermi surfaces split by the strong ASOC in the 
noncentrosymmetric compounds, which differs in origin 
from the usual T-independent van Vleck susceptibility 
stemming from the transition between the different or- 
bitals. A calculation assuming the simple spherical Fermi 
surface reveals that xf I'xS ~~ * 1 f° r H | [001] and 
xf c /x^ — ► 1/2 for H || [110] well below T c , irrespective of 
the pairing symmetry [4] , provided that the Rashba-type 
ASOC is sufficiently larger than the SC gap 15j. The ex- 
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perimental result is consistent in the case for H || [001], 
however, Xs°/Xs ~ 0.9 at T -> for H || [110] is not 
consistent with the prediction based on this model. Fu- 
jimoto pointed out theoretically that the X VV (T)/x s (T) 
in the SC state is intimately affected by electronic struc- 
tures and electron correlation effects [7|: In this case, the 
value of °/x^ a t low T limit corresponds to a fraction 
of x vv (T) contributing to Xs(T). Consequently, the re- 
sults of no change in Xa(T) for H || [001] across T c , but 
its slight decrease of Xs(T) for H || [110] below T c sug- 
gest the anisotropy of x VV (^)i dominated by the strong 
ASOC. This finding may be a characteristic feature of 
the spin susceptibility in the noncentrosymmetric super- 
conductor CeIrSi3 with the Rashba-type ASOC, which 
yields a large contribution of X VV (^) along the c-axis. 

This result also gives an indication for the large 
anisotropy of a Pauli limiting field (Hp) between H | 
[001] and H || [110] at T=0, which is roughly esti- 
mated through a comparison of a superconducting con- 
densation energy and a spin polarization energy given 
by Hp = k B T c /n B ^l - (x sc /x N ) 1- This fact gives a 
clue to understand why H C 2 is much larger along the c- 
axis than along the afo-plane in this compound. However, 
it should be noted that one possible reason why H C 2(0) 
in a narrow P range around P =2.6 GPa is exception- 
ally enhanced [13] may be attributed to strong AFM spin 
fluctuations at magnetic quantum criticality [23l ] or other 
unknown critical behaviors. Although the decrease in 
Xs(T) for H || [110] below T c evidences the Cooper-pair 
formation, it seems to be difficult to argue which types of 
Cooper-pairing states, that are theoretically predicted in 
noncentrosymmetric superconductors, take place: a spin- 
singlet state,a spin-triplet one, or an admixture of those 
states 0, H (! II H, [I] for the case that the strong 
ASOC dominates the spin susceptibility. 

Finally, we comment on the Knight-shift measure- 
ments in the HF superconductor CePtaSi without an in- 
version center along the c-axis Q: No change of ^Si- 
Knight shift for H parallel and perpendicular to the c- 
axis was observed in the samples either with T c =0.75 
K [ill or T c =0.46 K that is inherent to CePt 3 Si [H, [27, 
|28| . An anisotropy in the Knight shift was not resolved 
in CePtsSi, if any, because the lincwidth of 29 Si-NMR 
spectra was one order larger than in the single crystal of 
CeIrSi3 due to the coexistence of AFM order and the SC 
in CePt 3 Si at P=0. 

In conclusion, the extensive measurements of 29 Si 
Knight shift for the noncentrosymmetric CeIrSi3 under 
pressure have revealed that the spin component of the 
Knight shift if! 110 ' in the afc-plane slightly decreases, 
whereas kI° 01 ^ along the c-axis is unchanged across T c . 
This is considered because the spin susceptibility in the 
SC state is dominated by the van Vleck-like spin suscep- 
tibility in association with the strong ASOC caused by 
the absence of an inversion center as pointed out the- 



oretically. This is the first observation which exhibits 
an anisotropy of the spin susceptibility below T c in the 
noncentrosymmetric superconductor dominated by the 
strong Rashba-type spin-orbit interaction. We also com- 
ment that the Knight shift measurement does not give a 
direct evidence for the admixture of the Cooper-pairing 
state between even and odd parity in the case that the 
spin susceptibility is dominated by strong spin-orbit in- 
teraction. 
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